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FOREWORD 


This  report  describes  the  results  of  a  study  to  model  the  mass  distribution  of 
southern  Greenland  from  an  assortment  of  geophysical  data  for  the  purpose  of 
predicting  trends  in  the  geopotential.  The  work  was  performed  in  the  Space  and 
Surface  Systems  Division  of  the  Strategic  Systems  Department. 

This  report  has  been  reviewed  by  J.  L.  Sloop,  Head,  Space  and  Surface  Systems 
Division. 


Approved  by: 

R.  L.  SCHMIDT,  Head 
Strategic  Systems  Department 
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INTRODUCTION 


During  the  summer  seasons  of  1991  and  1992,  the  Naval  Research  Laboratory 
(NRL)  and  the  Naval  Oceanographic  Office  in  collaboration  with  researchers  from 
the  National  Oceanic  and  Atmospheric  Administration  and  the  National  Survey  and 
Cadastre  of  Denmark  will  perform  a  detailed  airborne  survey  of  Greenland  to  map 
gravity,  magnetics,  and  the  topography  of  the  entire  island.  During  the  first  summer 
campaign,  southern  Greenland,  below  70°  N,  will  be  surveyed  with  the  NRL  P-3 
Orion  aircraft  equipped  with  LaCoste-Romberg  and  Bell  gravimeters,  radar  and  laser 
altimeters,  Global  Positioning  System  (GPS)  receivers,  and  a  proton-precession 
magnetometer.  After  processing  of  the  GPS  data  to  remove  radial  accelerations  of 
the  aircraft,  it  is  anticipated  that  the  gravity  field  will  be  recovered  with  an  rms 
accuracy  of  2  mgal  for  half-degree  means.  Along-track  wavelength  resolution  should 
be  on  the  order  of  25-40  km.1 

Because  of  the  sparse  distribution  of  gravity  data  within  Greenland,  the  current 
mean  anomaly  field,  primarily  based  on  prediction^  fails  to  replicate  trends  in  the 
observed  geopotential  field.  To  better  forecast  trends,  a  project  was  initiated  at  the 
Naval  Surface  Warfare  Center  (NAVSWC)  to  collect  and  exploit  higher  resolution 
and  more  complete  geophysical  data  sources.  Such  improved  predictions  could  sup¬ 
port  survey  design  and  assist  in  the  interpretation  of  data  acquired  from  the  airborne 
survey.  The  result  of  this  effort  is  a  composite  topographic  data  base  that  supports 
the  prediction  of  intermediate  wavelengths  in  the  geopotential  field.  A  more  com¬ 
plete  model  of  the  geopotential  would  require  a  detailed  description  of  the  island's 
geology  which,  unfortunately,  is  well-known  only  in  the  coastal  ice-free  areas. 
However,  using  available  data  Haller3  has  postulated  a  general  description  of  the 
island's  geology. 

The  island  of  Greenland  is  part  of  the  Precambrian  shield  of  North  America. 
Two  Paleozoic  fold  belts  eroded  and  then  uplifted  by  more  recent  tectonic  activity  in 
the  coastal  areas  envelop  the  old  crystalline  block.  The  Paleozoic  plane  of  erosion  in 
the  center  of  the  island  may  be  underlain  by  two  northward-tilting  basement  blocks 
separated  by  a  distinct  depression  near  70°  N.  The  age  of  the  depression  may  be 
similar  to  the  age  of  the  coastal  tectonic  uplift.  More  recent  than  the  tectonic  uplift, 
isostatic  uplift  began  following  partial  removal  of  the  Pleistocene  ice  sheet.-1  The 
balance  of  that  ice  sheet,  which  averages  1.5  km  thick  over  most  of  the  island,3 
reaches  its  maximum  thickness  of  more  than  3100  m  near  72°  N.  Sizable  physio¬ 
graphic  features  such  as  these  make  significant  wavelength  contributions  to  the 
geopotential  model  because  their  gross  surfaces  are  well  defined  by  the  topographic- 
data  base. 


1 


NAVSWC  TR  91-409 


In  the  sections  that  follow,  the  data  sources  collected  during  the  project  are 
described.  The  later  sections  address  the  development  of  the  composite  topographic 
data  base,  geopotential  forecasting  using  the  inferred  mass  distribution,  and  model 
evaluations. 


GEOPHYSICAL  DATA  SOURCES 


Various  data  types  collected  for  the  Greenland  area  were  used  to  generate 
surface  elevations  of  the  ice  sheet  and  underlying  bedrock.  Other  data  sources  were 
used  to  evaluate  predictions  of  geopotential  trends  in  southern  Greenland.  Data 
sources  for  both  of  these  purposes  are  discussed. 


SATELLITE  RADAR  ALTIMETRY 

The  SEASAT  and  GEOSAT  satellite  radar  altimeters  were  designed  to  measure 
the  ocean  surfaces  to  recover  the  gravity  field  at  sea  and  mesoscale  oceanographic 
features.  New  techniques  developed  by  NASA  investigators^  enabled  the  use  of  these 
data  to  generate  elevations  of  non-ocean  surfaces  such  as  ice  sheets7  and  later  land 
areas. 8  Topographic  elevations  of  southern  Greenland  estimated  from  both  SEASAT 
and  GEOSAT  measurements  were  provided  by  NASA. 

Satellite  radar  altimeters  measure  range  between  the  orbit  and  a  reflecting 
surface  by  transmitting  a  signal  pulse  and  timing  the  delay  before  the  reflected  pulse 
is  received.  Successful  range  measurement  depends  on  accurate  onboard  prediction 
of  the  arrival  time  of  the  reflected  pulse.  If  the  reflected  pulse  arrives  within  the 
predicted  window,  tracking  is  maintained  and  onboard  range  estimates  are  made. 
These  ranges  are  telemetered  to  the  surface  along  with  the  total  reflected  waveform 
and  other  information.  When  the  reflected  pulse  arrives  outside  the  predicted  win¬ 
dow,  loss  of  track  occurs  and  no  range  estimates  are  formed.  Over  smooth  reflecting 
surfaces,  such  as  the  ocean,  the  reflected  pulse  assumes  a  relatively  simple  waveform 
and  range  estimates  made  onboard  are  accurate.  Over  rougher  surfaces,  such  as 
glacial  ice,  the  waveform  of  the  reflected  pulse  is  more  complex  and  gross  errors  in 
the  onboard  estimates  of  range  can  occur.  It  is  possible  to  make  better  estimates  of 
range  over  rougher  surfaces  by  reprocessing  the  telemetered  waveform.  This 
reprocessing  is  called  retracking M 

Retracking  a  telemetered  waveform  involves  fitting  it  iteratively  with  either  a 
five-parameter  or  nine-parameter  function,  depending  on  the  shape  of  the  waveform. 
Not  all  telemetered  waveforms  can  be  retracked  however,  and  those  that  cannot  be 
retracked  are  edited.  After  retracking  and  editing,  corrections  are  applied  to  the 
ranges  for  time  tag  biases,  signal  travel  time,  altimeter  antenna  offset,  ionospheric 
and  tropospheric  propagation  effects,  and  solid  Earth  tide  variations.  In  addition, 
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satellite  orbit  corrections,  primarily  in  the  radial  direction,  and  a  correction  for  a 
sloping  reflecting  surface1*1  are  also  applied.  These  corrected  ranges  are  used  to 
construct  elevations  of  the  reflecting  surface. 

Over  Greenland,  because  of  the  uneven  distribution  of  range  data  over  the  ice 
sheet,  regular  grids  of  elevations  have  been  developed.  11,12  Tangent  polar  stereo¬ 
graphic  projection  grids  have  been  generated  so  that  equal  spacing  between  grid 
points  occurs  on  the  projection  plane  tangent  at  the  North  Pole.  A  biquadratic  or 
bilinear  function  describing  the  surface  near  a  grid  point  was  fit  to  nearby  data.  Each 
altimeter  range  was  assigned  an  observation  standard  deviation  of  1.0  m.  After 
weighting  each  observation  by  the  inverse  of  distance  to  a  desired  grid  point,  a  least- 
squares  adjustment  was  made  so  that  data  closest  to  the  grid  point  influence  most  of 
its  estimated  elevation.  The  standard  deviation  of  a  grid  point  elevation  is  a  function 
of  how  well  the  elevation  represents  the  data. 

The  SEASAT  spacecraft  was  launched  in  late  June  1978  and  during  its  brief 
1 10-day  lifetime  collected  radar  altimetry  data  between  the  latitudes  72°  S  and  72°  N. 
Investigators  at  NASA  retracked  data  collected  over  southern  Greenland  and  used 
these  altimeter  ranges  to  generate  a  regular  grid  of  surface  elevations.1-  Approxi¬ 
mately  100,000  data  points  were  used  to  generate  a  tangent  polar  stereographic  grid 
with  spacing  between  grid  points  of  20  km.  Figure  1*  shows  the  geographic  pattern 
from  that  grid,  which  extends  over  the  ice  sheet,  and  exposed  land  of  some  coastal 
areas.  The  grid  elevations  were  referenced  to  the  GEM  10-B  geoid.  SEASAT-derived 
elevations  may  be  lower  on  average  than  elevations  derived  by  other  means  from 
between  2  to  3  mT3 

The  GEOSAT  spacecraft  was  launched  in  March  1985  and  provided  18  months 
of  high-resolution  radar  altimeter  data  between  latitudes  72°  S  and  72°  N  during  its 
primary  Geodetic  Mission.1*1  Eighteen  months  of  data  over  Greenland  were  re¬ 
tracked11  and  used  to  generate  two  grids.  Figure  2  shows  the  geographic  pattern 
from  the  polar  stereographic  grid  with  spacing  between  points  of  20  km.  Because 
GEOSAT  provided  as  much  as  10  times  the  amount  of  data  as  SEASAT  over  the  ice 
sheet  (personal  communication  with  Brenner,  1991),  NASA  also  generated  a  polar 
stereographic  grid  with  a  10-km  spacing  between  points.  The  elevations  in  both  grids 
were  referenced  to  the  GEM-10B  geoid. 

The  quality  and  quantity  of  non-ocean  measurements  from  the  two  altimeter 
missions  differ  because  of  hardware  and  software  modifications  to  GEOSAT.1 5  Over 
the  ice  sheet,  GEOSAT  provided  uenser  coverage  than  SEASAT.  Over  smaller  scale 
roughness  on  the  ice  s.ieet,  GEOSAT  tracked  better  because  of  modifications  to  its 
tracking  algorithm  However,  because  of  instrument  modifications,  once  GEOSAT 
lost  track  signal  reacquisition  was  not  as  quick  as  for  SEASAT.  For  this  reason,  over 
ice-free  rougher  terrain  where  tracking  loss  often  occurred,  SEASAT  provided  a 
better  description  of  the  topography. 


‘All  figures  arc  prescn ted  at  t  tic  end  of  t  tic  text 
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RADIO  ECHO  SOUNDINGS 

During  5  years  in  the  1970's,  1971-1972, 16  1974,17  1978, and  1979,11?  Danish 
investigators  collected  radio  echo  soundings  over  an  extensive  area  of  Greenland. 
The  data  were  collected  at  altitude  by  aircraft,  and  Figure  3  shows  the  tracks  of  data 
available  from  the  World  Data  Center- A  for  Glaciology  [Snow  and  Ice].  Nearly 
26,000  ice  surface  heights  and  nearly  18,000  underlying  bedrock  surface  heights 
were  available  for  southern  Greenland.  All  heights  are  relative  to  mean  sea  level  as 
measured  by  a  barometric  altimeter.  Figure  4  presents  an  example  of  800  nmi  of 
radio  echo  soundings  collected  in  1974  between  68°  and  74°  N.  In  this  figure  the  ice 
surface  occasionally  reaches  more  than  2800  rn  above  the  bedrock  surface. 

The  principle  of  a  radio-echo  sounder  is  described  by  Evans  and  Smith. Dis¬ 
tances  between  the  instrument  and  reflecting  surface  are  measured  by  transmitting 
a  radar  pulse  and  timing  the  delay  before  the  reflected  pulse  is  received.  Measuring 
distances  to  the  ice  and  ice-covered  bedrock  surfaces  requires  two  signals  of  different 
frequencies.  The  higher  frequency,  300  MHz,  is  reflected  by  the  ice  surface  while  the 
lower  frequency,  60  MHz,  penetrates  the  ice  and  is  reflected  by  the  bedrock.  During 
data  collection,  oscilloscope  images  of  the  reflected  signals  are  recorded  on  35mm  film 
for  later  analysis. 

The  images  on  the  film  give  a  measure  of  the  time  elapsed  between  signal 
transmission  and  echo  return  subject  to  propagation  variations  and  systematic 
errors.  For  all  bedrock  soundings,  the  electromagnetic  wave  velocity  of  169  m/ps  in 
solid  ice  with  a  density  of  920  kg/m3  was  assumed.  Though  wave  velocity  varies  with 
ice  density,  no  corrections  were  made  because  the  density  profile  of  the  ice  sheet  is 
generally  unknown  because  of  the  limited  distribution  of  ice  cores. 31, 32  In  addition, 
the  broad  radiation  pattern  of  the  antenna  may  result  in  reflections  from  surfaces  not 
directly  below  the  antenna.  Landforms,  such  as  mountains,  or  inhomogeneities  in 
the  ice  sheet,  such  as  ice  lenses  in  the  firn,  far  away  from  the  aircraft,  may  appear  on 
the  sounder.  Though  this  phenomenon  may  add  noise  to  measurements,  no  correc¬ 
tions  were  made.  Ranges  were  corrected  for  the  cable  length  between  the  transmitter 
and  the  antenna  because  transmission  time  was  measured  from  the  transmitter.  The 
resulting  range  accuracy  of  the  radar  is  10  m  or  1  percent  of  the  measured  depth, 
whichever  is  greater. 30  After  computing  ranges,  the  positions  of  the  aircraft  were 
used  to  determine  mean  sea  level  elevations. 

Associated  with  each  measurement  are  coordinates  determined  by  an  Inertial 
Navigation  System  (INS)  onboard  the  aircraft.  Coordinates  determined  by  the  three- 
axis  INS  have  a  horizontal  accuracy  of  1-5  km. 33  Inaccuracies  in  aircraft  altitude 
were  large  and  required  correction.  Locations  defined  by  the  Navy  Navigation 
Satellite  System  (NNSS)  and  other  sources  provided  reference  points  to  calibrate 
track  altitudes.  Crossovers  with  calibrated  tracks  were  used  to  correct  those  tracks 
not  passing  over  reference  points.  Additional  comments  on  the  vertical  accuracy  of 
the  aircraft  positions  are  provided  below. 
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GEOLOGICAL  MAP 

A  tectonic/geological  map  of  Greenland, ^4  developed  by  the  Geological  Survey  ol 
Greenland,  was  obtained  through  the  Library  of  Congress.  Of  interest  on  this  map 
are  contours  of  topographic  heights  of  the  ice  sheet  and  exposed  land  near  the  coast. 
Also  contoured  is  the  surface  underlying  the  ice.  The  sources  of  this  topographic 
information  were  maps  from  the  Geodaetisk  Institut,  Kobenhaven,  and  World 
Aeronautical  Charts  published  by  the  U.  S.  Air  Force. 


GRAVITY  ANOMALIES 

Gravity  measurements  collected  near  the  Greenland  coast  by  surface  ship  and 
overland  by  various  means  were  compiled  by  NRL  and  the  Defense  Mapping  Agency 
Aerospace  Center  (DMAAC).  Many  sources  contributed  to  these  data,  but  the  pri¬ 
mary  sources  were  from  Greenland,  Denmark,  the  United  States,  Canada,  and 
Sweden.  The  distribution  of  data  compiled  by  NRL  is  shown  in  Figure  5  and  by 
DMAAC  in  Figure  6.  As  can  be  seen  from  these  two  Figures,  significant  overlap  exists 
between  the  overland  data  held  by  these  two  organizations.  Further,  it  is  evident 
that  an  abundance  of  offshore  data  exists.  These  data  are  primarily  sea  surface 
measurements  with  an  occasional  sounding.  Data  on  land  are  concentrated  sporad¬ 
ically  near  the  coast  with  four  notable  ice-traversing  tracks.  In  both  data  sets,  one 
track  appears  near  75°  N  with  two  more,  nearly  overlapping  tracks  around  7^°  N.  In 
the  DMAAC  data  set  a  short  45-km  track  near  67°  N  also  appears.  The  overland 
track  data  were  collected  by  an  aircraft  landing  on  the  ice.  The  data  consist  of  the 
observation  station  coordinates  and  elevation,  observed  gravity,  and  free-air  and/or 
Bouguer  anomaly. 

The  gravity  anomaly  describes  perturbations  of  the  field  caused  by  variations  in 
mass  distribution.  It  is  the  difference  between  observed  gravity  reduced  to  the  geoid 
and  a  theoretical  value  on  the  reference  ellipsoid.  The  international  gravity 
formula 25  predicts  the  value  of  gravity  on  the  ellipsoid  as  a  function  of  latitude. 
Various  corrections  are  introduced  to  reduce  the  observed  value  to  the  geoid.  The 
free-air  correction  adds  to  the  observed  value  the  effect  of  the  elevation  of  the 
observation  station.  The  Bouguer  correction  is  introduced  to  account  for  the  masses 
between  the  station  and  sea  level.  The  purpose  for  the  gravity  reduction  defines  the 
type  of  corrections  that  are  applied  to  the  observed  gravity  value.  In  the  analysis 
presented  below,  free-air  anomalies  were  used  because  of  their  strong  correlation 
with  topography. 
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METHODS 


Using  these  data  sets,  surface  elevation  grids  of  the  ice  and  bedrock  topography 
of  southern  Greenland  were  developed.  The  grids  were  used  to  model  the  geopoten¬ 
tial  of  the  same  area. 


SURFACE  ELEVATION  GRIDS 

Independently  no  single  data  set  sufficiently  describes  either  the  ice  or  bedrock 
surface.  For  this  reason,  composite  data  sets  for  the  ice  surface  and  bedrock  surface 
were  developed.  This  involved  editing  data  sets  based  on  their  relative  strengths  and 
merging  them.  After  the  composite  data  sets  were  developed,  they  were  interpolated 
into  rectangular  grids. 

The  topographic  elevations  of  Greenland  south  of  72°  N,  including  the  inland  ice 
sheet  and  coastal  ice-free  areas,  will  be  referenced  as  the  composite  ice  surface.  The 
composite  ice  surface  was  developed  from  three  data  sources.  The  GEOSAT  20-km 
elevation  grid  provided  data  for  the  interior  ice  sheet,  which  were  consistent  with 
large-scale  features  found  on  the  geological  map.  Further,  comparisons  of  these  data 
with  the  SEASAT  20-km  grid  elevations26  indicated  a  significant  average  height 
increase,  which  is  consistent  with  a  growing  ice  sheet.  However,  in  coastal  ice-free 
land  areas  defined  by  Bindschadler's  maps, 13  the  GEOSAT  20-km  grid  contained 
sporadic  negative  elevations.  These  were  inconsistent  with  elevation  contours  on  the 
geological  map  and  different  from  the  SEASAT  20-km  grid  elevations*  by  as  much  as 
1.99  km.  For  these  reasons,  the  GEOSAT  negative  elevations  were  manually  edited 
and  replaced  with  SEASAT  or  other  elevation  data  (Figure  7).  SEASAT  elevations 
were  substituted  when  consistent  with  the  geological  map  and  when  SEASAT  grid 
point  uncertainties  were  significantly  smaller  than  those  from  GEOSAT.  In  loca¬ 
tions  where  neither  GEOSAT  nor  SEASAT  provided  reliable  data,  elevations  were 
estimated  using  averages  of  nearby  elevations  from  the  NRL  gravity  data  base.  The 
remaining  negative  heights,  principally  near  the  coast,  were  substituted  with  values 
interpolated  between  the  ocean  surface  and  valid  data  base  grid  points  above  sea 
level. 


The  elevations  of  ice-covered  bedrock  and  coastal  ice-free  areas  south  of  72°  N 
will  be  referenced  as  the  composite  bedrock  surface.  The  composite  bedrock  surface 
was  developed  from  the  Danish  radio  echo  soundings  of  the  bedrock  and  other  data. 
In  the  coastal  ice-free  areas  where  gaps  exist  in  the  radio  echo  soundings,  the  Danish 
data  were  supplemented  with  the  manually  edited  elevations  from  the  composite  ice 
surface  data  base,  their  locations  shown  in  Figure  7.  These  data  were  used  to 
generate  the  composite  bedrock  surface  elevation  grid. 

*  In  I  he  coastal  areas  where  terrain  is  rougher  than  over  the  inland  ice,  SEASAT  generally  described 
the  topography  better  t  han  GEOSAT  (see  the  Satellite  Radar  Altimetry  section). 
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In  addition  to  the  composite  ice  surface  elevation  grid  developed  primarily  with 
satellite  altimetry,  a  second  grid  of  that  surface  was  generated.  Using  the  Danish 
soundings  of  the  ice  surface  and  the  same  manually  edited  elevations  supplementing 
the  bedrock  elevations,  the  Danish  ice  surface  elevation  grid  was  generated. 

All  three  surface  elevation  grids  were  developed  from  their  source  data  in  a 
similar  manner.  Each  respective  data  set  was  augmented  with  zeros  representing 
ocean  surface  elevations.  These  zeros  extended  outward  from  the  shoreline,  making 
the  dimensions  of  the  rectangular  grid  area  38.4  longitude  degrees  by  12.0  latitude 
degrees.  Using  these  augmented  data  sets,  grids  with  a  spacing  between  grid  points 
of  0.2  deg  in  latitude  and  longitude  were  developed. 

A  technique  developed  by  Braile-7  was  used  to  interpolate  the  data  onto  a  regu¬ 
lar  rectangular  grid.  This  technique  is  based  on  fitting  a  polynomial  to  a  small  region 
of  the  grid  area.  The  fitting  region  contains  several  grid  points  and  is  moved  succes¬ 
sively  as  a  template  over  the  entire  grid  area.  The  characteristics  of  the  gridding 
operation  can  be  controlled  by  the  choices  of  the  order  of  the  polynomial,  the  size  of 
the  fitting  region,  and  the  number  of  grid  points  to  interpolate  in  each  fit  region.  The 
best-fitting,  in  a  least-squares  sense,  surface  is  determined  and  from  it  grid  points  are 
interpolated. 

The  order  of  the  polynomial  used  is  dependent  on  the  number  of  data  points  in 
the  fitting  region.  The  chosen  order  of  the  polynomial  is  reduced  dynamically  until 
the  number  of  points  available  are  sufficient  for  solution.  This  provides  for  optimum 
fitting  in  areas  of  high  and  low  data  density  and  minimizes  edge  effects.  The  grids  of 
all  three  surfaces  were  made  using  a  third-order  polynomial.  The  fit  region  template 
for  the  ice  surface  grid  was  2  deg  square  while  those  of  the  bedrock  and  Danish  ice 
surface  grids  were  3  deg  square.  Larger  fitting  regions  were  used  for  the  latter  two 
surfaces  because  of  track-spacing  irregularity  and  gaps  in  data  distribution.  Sixteen 
grid  points  were  interpolated  per  fit  region  r  :  all  three  surfaces. 

The  first  two  of  these  grids  were  subsequently  used  to  generate  a  model  of  the 
geopotential  trends  in  southern  Greenland.  The  third  grid,  the  Danish  ice  surface, 
was  used  to  evaluate  the  elevation  accuracy  of  the  bedrock  surface  where  defined  by 
the  radio  echo  sounding  data. 


GEOPOTENTIAL  MODELING 

The  composite  topographic  data  base  of  southern  Greenland  consists  of  mean  sea 
level  elevations  of  the  ice  surface  and  bedrock  surface  in  grids  with  a  spacing  between 
points  of  0.2  deg.  This  data  base  provides,  assuming  a  density  structure,  the  localized 
mass  distribution  and  thus  insight  into  geopotential  components.  As  constructed 
from  available  source  data,  the  elevation  grids  cover  the  area  between  60.5°  and  72.3° 
N  latitude  and  between  20.8°  and  59.0°  W  longitude.  There  are  11520  data  base 
elements  each  having  two  elevation  values.  Elements  located  offshore  had  zero 
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elevations  before  the  gridding  procedure;  however,  because  of  edge  effects  of  the 
surface  fit,  small  nonzero  elevations  exist  in  some  locations  near  shore.  For  inland 
elements,  the  lower  elevations  define  either  the  exposed  bedrock  surface  or  the  inter¬ 
face  between  the  bedrock  and  overlying  ice.  The  upper  elevations,  when  different, 
provide  the  surface  topography  of  the  Greenland  ice  sheet.  Thus,  the  composite 
topographic  data  base  defines  the  mass  distribution  for  southern  Greenland  as  two 
sets  of  rectangular  prisms  with  common  element  bases,  see  Figure  8. 

The  geopotential  effect  of  this  irregular  mass  distribution  can  then  be  calcu¬ 
lated  at  discrete  points.  Assuming  a  constant  density  for  the  material  within  each 
prism  element,  the  data  base  becomes  a  finite  element  representation  for  the  mass 
distribution.  Taking  2.67  g/cm^  as  the  density  of  the  bedrock  layer  and  1.00  g/cm'i 4 
for  the  ice  layer,  two  integrations  over  the  mass  layers  are  required  to  calculate  the 
effect  of  the  masses  at  a  given  point. 

The  numerical  integration  is  accomplished  using  the  methods  described  in 
Forsberg,^  which  provides  closed-form  expressions  for  the  contribution  of  individual 
elements  (rectangular  prisms)  to  various  gravimetric  quantities.  At  each  computa¬ 
tion  point,  these  formulas  are  evaluated  then  summed  over  all  elements  of  the  com¬ 
posite  data  base,  implying  up  to  23040  computations  for  each  point,  since  two  distinct 
density  layers  are  present.  In  practice,  templates  of  desired  radii  may  be  used  to 
reduce  the  complexity  of  the  computations.  An  inner  template  defines  the  region 
where  exact  formulas  are  used.  A  larger  template  zone  is  defined  in  which  approxi¬ 
mations  are  made  to  increase  computation  speed.  Prisms  external  to  the  second 
template  radius  are  processed  with  still  further  approximations.  The  template  radii 
for  the  computations  performed  in  this  study  were  50  and  1000  km. 

The  result  of  this  integration  is  the  components  of  various  gravimetric 
quantities  due  to  the  anomalous  mass  distribution  above  mean  sea  level.  These 
components  are  wavelength  limited  by  the  resolution  of  the  data  base  and  represent 
only  part  of  the  amplitude  of  the  full  geopotential  field.  The  results  of  modeling  only 
the  anomalous  mass  distribution,  if  accurate,  should  correlate  with  intermediate 
wavelengths  of  the  full  geopotential  field  and  provide  insight  into  geopotential  trends 
in  the  region.  Geopotential  forecasts  based  on  this  model  are  compared  with  free-air 
anomaly  profiles  below. 


1  1  00  (’/cm3  was  used  m  all  calculal  ions  involving  the  ice  sheet.  Although  the  density  of  ice  at  0°  C  is 

0.92  gem3,  the  temperature  and  pressure  distribution  throughout  the  Greenland  ice  sheet  is 
complicated  and  unknown.  Gi  ven  the  resolution  and  accuracy  of  the  data  base,  the  assumed  density 
was  considered  sufficient  lor  I  he  purpose  of  this  analysis. 


8 


NAVSWC  TR  91-409 


RESULTS 


COMPOSITE  TOPOGRAPHIC  DATA  BASE 

The  contoured  ice  surface  presented  in  Figure  9,  and  its  three-dimensional 
perspective  in  Figure  10,  show  features  similar  to  those  generated  by  NASA  inves¬ 
tigators.!  M3  The  highest  elevations,  over  3100  m,  occur  near  the  northern  extent  of 
the  area.  An  ice  ridge  strikes  NNE  through  the  center  of  the  ice  sheet  flanked  by 
steep  gradients  on  the  east  coast  and  more  ge  ’e  gradients  on  the  west. 

The  contoured  bedrock  surface  presented  in  Figure  11,  and  its  three- 
dimensional  perspective  in  Figure  12,  show  physiographic  features  consistent  with, 
but  in  much  more  detail  than,  the  geological  map.  Apparent  on  the  east  coast  are  the 
steep  gradients  of  mountain  ranges  formed  by  uplifted  Paleozoic  folds.  In  the  center 
of  the  island,  the  southern  northward-tilting  basement  block  may  be  evidenced  by  the 
slope,  in  that  direction.  The  depression,  possibly  separating  the  two  tilted  basement 
blocks,  can  be  seen  between  68°  N  and  70°  N. 

The  contoured  Danish  ice  surface  presented  in  Figure  13  has  features  similar  to 
the  ice  surface  in  Figure  9.  The  mean  elevation  difference*  between  these  two  sur¬ 
faces  is  -12  m  with  a  sigma  of  170  m.  The  sign  of  this  result  indicates  that  the  com¬ 
posite  ice  surface  is  on  the  average  slightly  higher  than  the  Danish  ice  surface. 
Figure  14  gives  the  locations  of  differences  whose  absolute  values  were  greater  than 
100  m.  These  differences  are  primarily  located  in  the  coastal  regions,  where  gaps  in 
the  Danish  soundings  exist,  and  near  the  northern  limit  of  the  map  area  where  grid 
edge  effects  may  introduce  variability.  Although  differences  between  these  two 
surfaces  exist,  the  Danish  ice  surface  is  generally  consistent  with  the  composite  ice 
surface.  This  is  evidence  that  the  vertical  positions  of  the  aircraft  collecting  these 
data  were  well  determined,  and  the  Danish  ice  surface  soundings  are  an  accurate 
data  set.  Further,  the  Danish  bedrock  surface  soundings,  collected  simultaneously 
with  the  ice  surface  soundings  and  reduced  in  a  similar  manner,  have  acceptable 
vertical  control.  This  suggests  that  the  composite  bedrock  surface  is  generally  an 
accurate  description  of  the  physical  bedrock  surface  where  data  are  available. 


GRAVITY  DISTURBANCE  AND  VERTICAL  DEFLECTION 
TRENDS  IN  SOUTHERN  GREENLAND 

With  assumed  values  for  the  density  of  the  ice  and  bedrock  layers,  the 
topographic  data  base  was  used  to  approximate  the  mass  distribution  of  southern 


*  On t *  (i if ierem es  whose  absolute  values  were  greater  than  0  05  m  were  used  to  determine  these 
st  at  1st  ies,  thus  ell  minat  ocean  },rrid  points  with  null  e  lev  a  I  ions 
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Greenland.  The  geopotential  effect  of  this  irregular  mass  distribution  was  then 
calculated  for  discrete  points  lying  external  to  the  mass.  The  quantities  determined 
through  these  computations  were  gravity  disturbance  and  vertical  deflection. 

Since  the  survey  of  Greenland  is  planned  to  he  flown  at  constant  altitude,  geo¬ 
potential  computations  for  southern  Greenland  were  made  on  a  constant  altitude 
surface  2000  ft  above  the  highest  grid  point  in  the  data  base  (3195  m).  The  computed 
values  of  gravity  disturbance  and  vertical  deflection,  at  a  grid  spacing  equivalent  to 
that  of  the  data  base,  were  contoured.  The  results  are  provided  in  Figure  15  for 
gravity  disturbance  and  in  Figures  16  and  17  for  the  two  components  of  vertical 
deflection.  The  gravity  disturbances  at  altitude  provide  insight  into  the  spatial 
wavelengths  of  the  Field  and  thus  can  contribute  to  the  airborne  survey  design.  For 
example,  in  the  central  longitude  zones  of  the  island,  where  the  ice  sheet  is  the 
dominate  feature,  the  geopotential  field  is  more  slowly  varying.  To  the  resolution 
and  accuracy  of  the  data  base,  there  are  east-west  variations  as  small  as  70  mgal 
over  distances  as  great  as  15  deg.  Contrasted  to  this  are  the  coastal  zones,  which 
demonstrate  variations  in  gravity  disturbance  of  up  to  200  mgal  over  2.5-deg  arc 
lengths,  with  even  higher  gradients  in  some  regions.  These  high-amplitude  spatial 
variations  are  associated  with  the  large  variations  in  the  bedrock  topography  in 
many  coastal  areas  (see  for  instance  Figure  11).  These  larger  gradients  are  also 
found  in  the  east-west  component  of  vertical  deflection  along  the  eastern  coastal 
regions  as  shown  in  Figure  17. 

A  three-dimensional  perspective  of  the  disturbances  at  altitude  is  given  by 
Figure  18.  Notice  that  the  disturbance  field  is  highly  correlated  with  the  bedrock 
topography  (Figure  12)  as  expected. 

In  terms  of  survey  design  for  geopotential  recovery,  these  results  suggest  that 
denser  spaced  tracks  are  required  in  an  east-west  direction  in  many  coastal  regions 
with  greater  track  spacing  being  adequate  over  most  of  the  central  ice  sheet.  Without 
the  use  of  a  track-spacing  factor  that  is  geographically  adjusted,  it  is  likely  that  the 
mean  gravity  anomalies  recovered  from  the  airborne  data  will  demonstrate 
significant  biases  in  these  high  gradient  provinces. 

In  the  section  that,  follows,  examples  of  gravity  disturbances  on  the  surface 
computed  from  the  topographic  model  will  be  compared  with  measured  free-air 
anomalies. 


COMPARISON  OF  GEOPOTENTIAL  FORECASTS  WITH 
FREE-AIR  ANOMALY  DATA 

Linear  Profiles 


The  gravity  disturbance  on  the  physical  surface  P'  is  related  to  the  free-air 
gravity  anomaly  to  first  order  by  the  following  equation: 
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where  P  is  the  corresponding  point  on  the  geoid;  H  is  orthometric  height;  N  is  the 
undulation  at  P  and  y  is  normal  gravity.  Examining  the  contribution  of  each  term  in 
the  right-hand  side  of  Equation  1,  it  is  evident  that  the  gravity  disturbance  will  be 
spatially  correlated  with  free-air  anomaly. 

In  the  comparisons  that  follow,  that  part  of  the  total  gravity  disturbance  of 
Equation  1  due  to  the  anomalous  mass  distribution  above  mean  sea  level  is  computed 
from  the  composite  topographic  data  base  at  the  surface  locations  P'  where  point 
anomaly  data  have  been  measured.  Since  this  component  of  the  total  disturbance  is 
due  to  the  anomalous  mass,  it  is  expected  to  also  demonstrate  a  strong  correlation 
with  topography  and  thus  with  the  free-air  anomaly,  again  to  the  resolution  and 
accuracy  of  the  data  base. 

Figure  19  shows  the  location  of  a  primarily  east-west  anomaly  profile  across 
southern  Greenland  found  in  the  NRL  holdings.  The  profile  has  a  longitude  extent  of 
over  16  deg.  Using  the  data  base  model,  gravity  disturbances  were  computed  at  the 
locations  of  these  surface  grav  .ty  stations.  The  location  of  each  point  was  extracted 
from  the  NRL  gravity  anomaly  file  along  with  the  point  elevation.  Figure  20 
provides  the  free-air  anomaly  values  along  the  profile  and  the  gravity  disturbances 
simulated  from  the  data  base  model.  The  disturbances  here  have  been  translated  by 
the  subtraction  of  93.1  mgal  to  bring  the  two  scales  into  general  agreement.  A 
comparison  of  these  profiles  shows  that  the  overall  trend  in  the  actual  anomaly  data 
appears  in  the  simulated  results.  Several  of  the  shorter  wavelength  variations  are 
common  to  both  profiles  although  the  variation  in  magnitude  of  the  simulated  data  is 
much  less. 

Figure  21  gives  the  location  of  another  free-air  anomaly  profile  found  in  the 
data  provided  by  DMAAC.  This  profile  covers  approximately  2.5  deg  of  longitude 
near  the  eastern  coast.  The  disturbances  at  these  locations  predicted  by  the  topo¬ 
graphic  model,  translated  by  -79.7  mgal,  are  found  in  Figure  22  along  with  the 
anomaly  data.  Again  the  overall  decreasing  trend  in  the  free-air  data  is  present  in 
the  disturbance  profile. 

Figure  23  shows  a  companion  profile  from  the  NRL  data  set  to  the  profile  in 
Figure  19,  with  approximately  15  deg  extent  in  longitude.  Figure  24  provides  the 
fiee-air  anomaly  and  computed  disturbance  profiles  after  removal  of  83.0  mgal.  In 
this  example,  the  general  trend  over  the  extent  of  the  track  is  present  in  the  dis¬ 
turbance  results;  however,  its  slope  is  significantly  greater  than  in  the  anomaly  data 
Shorter  wavelengths  seen  in  the  anomaly  data  are  generally  missing  from  the 
disturbance. 
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Two  points  need  to  be  made  regarding  these  results.  First,  using  the  topo¬ 
graphic  data  base  and  constant  layer  density,  intermediate  wavelength  trends  in  the 
observed  anomaly  data  can  generally  but  not  always  be  predicted.  Because  of  the 
resolution  and  accuracy  of  the  data  base  and  density  assumptions,  short  wavelengths 
in  the  Field  are  generally  not  represented  by  the  model.  The  shortest  wavelengths 
that  can  be  represented  are  about  0.4  deg.  Thus,  the  computed  disturbance  Field  is 
generally  smoother  than  what  will  be  acquired  by  the  airborne  survey.  Second,  since 
the  disturbances  are  based  solely  on  the  modeled  representation  of  anomalous  mass 
above  mean  sea  level,  all  numerical  results  will  be  biased  due  to  the  geopotential  of 
unmodeled  mass  below  mean  sea  level. 


Area  FroFiles 


To  further  evaluate  the  model,  three-dimensional  perspectives  of  observed  free- 
air  anomalies  and  computed  gravity  disturbances  were  generated.  The  observations, 
collected  on  the  surface  in  several  regions  near  the  coast,  are  irregularly  distributed 
and  contain  gaps.  The  gravity  disturbances  were  computed  at  the  same  locations 
where  observations  were  collected.  Before  perspectives  were  generated,  regular  grids 
were  developed  during  which  values  were  interpolated  in  areas  void  of  data. 
Therefore,  the  features  in  these  areas,  shaded  in  the  Figures,  are  probably  inaccurate 
and  should  be  ignored.  The  coastal  regions  presented  often  coincide  with  areas  where 
data  base  elevations  are  weakest  and  where  extensive  data  editing  was  necessary. 
For  these  reasons,  some  disagreement  is  expected  in  the  comparisons  between 
perspectives  of  the  computed  gravity  disturbances  and  those  produced  from  measured 
gravity  anomalies. 

The  First  area  considered,  known  as  Gasefjord,  is  shown  in  Figure  25.  Twenty- 
Five  measured  free-air  anomalies  in  the  area  were  found  in  the  gravity  data  provided 
by  NRL.  An  enlargement  of  the  area  showing  the  locations  of  the  free-air  anomalies 
is  given  in  Figure  26.  Figure  27  shows  the  three-dimensional  perspective  of  these 
anomalies  and  the  computed  gravity  disturbances.  In  the  rectangular  area 
presented,  there  is  no  occurrence  of  measured  anomalies  (see  Figure  26)  in  the 
southv'est  corner  so  the  corresponding  area  is  shaded  in  Figure  27.  A  comparison  of 
the  light  surfaces  shows  some  common  trends.  For  instance,  each  surface  shows  a 
dominant  crest  extending  west  from  near  69.5°  N.  On  the  southwest  side  near  27.°  W, 
a  smaller  convex  feature  appears  on  both  surfaces. 

Figure  28  provides  the  value  of  the  free-air  anomaly  and  disturbance  at  each  of 
the  25  points  of  Figure  26  numbered  consecutively  from  west  to  east.  Of  the  24  con¬ 
secutive  pairs  of  points  in  this  ordering  scheme,  the  direction  (increase  or  decrease)  of 
change  in  the  anomaly  was  predicted  using  the  model  disturbances  in  22  cases. 
Success  in  predicting  the  magnitude  of  the  change  varied. 

A  second  example  is  provided  by  the  anomaly  data  available  in  Holsteinsborg, 
Figures  29  and  30.  The  three-dimensional  perspective  of  the  anomaly  data  and 
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modeled  disturbance  appear  in  Figure  31.  Again  trends  in  the  anomaly  data  are 
predictable  using  the  modeled  disturbances.  For  example,  a  prominent  eastward 
slope  appears  in  the  anomalies  and  the  disturbances. 

A  third  example  of  this  type  of  comparison  is  Jameson  Land,  Figures  32,  33, 
and  34. 


CONCLUSIONS 


Regular  grids  of  the  ice  surface  and  bedrock  surface  of  southern  Greenland  with 
assumed  density  structures  provide  models  of  the  localized  mass  distribution  and 
insight  into  intermediate  wavelengths  of  the  geopotential.  The  geopotential  model  is 
limited  by  the  resolution  and  accuracy  of  the  composite  topographic  data  base  and  the 
assumption  of  density  homogeneity  in  the  ice  and  bedrock  layers.  Track  and  regional 
comparisons  between  the  modeled  geopotential  and  free-air  anomaly  data  collected 
on  the  island  generally  show  similarities  in  their  intermediate  wavelengths.  Though 
short  wavelengths  of  the  anomalous  field  are  not  represented  in  the  model,  it  can 
support  survey  design  and  interpretation  of  new  survey  data  from  the  region. 
Further  comparisons  of  the  model  with  data  acquired  in  1991  are  planned. 
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FIGURE  3.  TRACK  DISTRIBUTION  OF  DANISH  RADIO  ECHO  SOUNDINGS 
SHOWN  ARE  5  YEARS  OF  DATA. 


FIGURE  4.  EXAMPLE  OF  DANISH  RADIO  ECHO  SOUNDINGS  PROFILE 
COLLECTED  IN  1974  BETWEEN  68°  N  AND  74°  N. 
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FIGURE  7.  LOCATION  a  OF  MAjN  UALLY  EDITED  DATA 


FIGURE  8.  DUAL  PRISM  MODELS  WITH  COMMON  BASES 
FOR  ICE  AND  BEDROCK  MASS  DISTRIBUTIONS 


19 


NAVSWCTR  91-409 


20 


NAVSWCTR  91-409 


21 


FIGURE  11.  CONTOUR  OF  BEDROCK  ELEVATIONS  FOR  FIGURE  12.  3D  PERSPECTIVE  OF  BEDROCK  ELEVATIONS 

SOUTHERN  GREENLAND.  CONTOUR  INTERVAL  IS  100  M.  OF  SOUTHERN  GREENLAND.  VIEW  IS  FROM  SSW 

VERTICAL  SCALE  EXAGGERATED 
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FIGURE  19.  NRL  GRAVITY  ANOMALY  PROFILE 
NEAR  70.4°  N 


FIGURE  20,  ACTUAL  ANOMALY  VS  SIMULATED 
DISTURBANCE  ALONG  NRL  PROFILE 
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FIGURE  21.  DMAAC  GRAVITY  ANOMALY  PROFILE 
NEAR  66  4°  N 


FIGURE  22  ACTUAL  ANOMALY  VS  SIMULATED 
DISTURBANCE  ALONG  DMAAC  PROFILE 
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FIGURE  23.  NRL  GRAVITY  ANOMALY  PROFILE 
NEAR  70.5°  N 
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FIGURE  24  ACTUAL  ANOMALY  VS  SIMULATED 
DISTURBANCE  ALONG  NRL  PROFILE 
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NRI  FRKK-A1R  ANOMA1.IKS  GRAVITY  D1STURBANCK  BV  TOPOGRAPHY 


FIGUILc  ..  GASEFJORD  3»  PERSPECTIVE  OF  FREE-AIR  AND 
SIMUL/  '.F  GRAVITY.  VIEW  IS  FROM  SSE  SHADED  AREAS  LACK  DATA 
AND  ARE  PROBABLY  INACCURATE.  VERTICAL  SCALE  I  =  15  MGAL. 


POINT  NUMBER 

FIGURE  28.  MODEL  COMPARISON  IN  GASEFJORD  AREA 
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FIGURE  29.  HOLSTEINSBORG  LOCATION  OF  FIGURE  30.  HOLSTEINSBORG  DISTRIBUTION  OF 

NRL  FREE-AIR  ANOMALIES  NRL  FREE-AIR  ANOMALIES 
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FIGURE  31.  HOLSTEINSBORG  3»  PERSPECTIVE  OF  FREE-AIR  AND 
SIMULATED  GRAVITY.  VIEW  IS  FROM  ESE.  SHADED  AREAS  LACK  DATA 
AND  ARE  PROBABLY  INACCURATE.  VERTICAL  SCALE  I  =  7  MGAL. 


FIGURE  32  JAMESON  LAND  LOCATION  OF  NRL  FREE-AIR  ANOMALIES 
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FIGURE  33.  JAMESON  LAND  DISTRIBUTION  OF 
NRL  FREE-AIR  ANOMALIES 


NRL  FREE-AIR  ANOMALIES  GRAVITY  DISTURBANCE  BY  TOPOGRAPHY 


FIGURE  34  JAMESON  LAND  3*>  PERSPECTIVE  OF  FREE-AIR 
AND  SIMULATED  GRAVITY  VIEW  IS  FROM  WEST. 
VERTICAL  SCALE  1  =  15MGAL. 


32 


NAVSWCTR  90-409 


DISTRIBUTION 


Copies 

ATTN  K  BURKE  (PPG)  1 

W  WOODEN  (PRAi  1 

DIRECTOR 

DEFENSE  MAPPING  AGENCY 

HEADQUARTERS 

BLDG  56 

US  NAVAL  OBSERVATORY 
WASHINGTON  DC  20305-3060 


ATTN  L  DOTSON  iDSG) 
DIRECTOR 

DEFENSE  MAPPING  AGENCY 
AEROSPACE  CENTER 
3200  S  SECOND  ST 
ST  LOUIS  MO  63118-3399 


ATTN  R  ZIEGLER  (SGS)  5 

M  KUMAR  tSGS)  1 

DIRECTOR 

DEFENSE  MAPPING  AGENCY 
SYSTEMS  CENTER 
STE  800 

8301  GREENSBORO  DR 
MCLEAN  VA  22102 


ATTN  J  BROZENA  5 

M  PETERS  1 

P  VOGT  1 

NAVAL  RESEARCH  LABORATORY 
MAIL  CODE  51 10 
4555  OVERLOOK  AVE  SW 
WASHINGTON  DC  20375 


ATTN  S  BARKER  5 

M  CHALONA  1 

COMMANDING  OFFICER 
NAVAL  OCEANOGRAPHIC  OFFICE 
CODE  GG 

STENNIS  SPACE  CENTER  MS  39522-5001 


Copies 

ATTN  G  MADER  1 

N  GSN  OS'  NO  A  A  N/CG114 
6001  EXECUTIVE  BLVD 
ROCKVILLE  MD  20852 


ATTN  R  RAPP  1 

DEPT  OF  GEODETIC  SCIENCE 

AND  SURVEYING 

THE  OHIO  STATE  UNIVERSITY 

1958  NEIL  AVE 

COLUMBUS  OH  43210 


ATTN  CODE  671 

R  B1NDSCHADLER  1 

H  ZWALLY  1 

CODE  626 

J  MARSH  1 

CODE  920 

H  FREY  1 


NATIONAL  AERONAUTICS  AND 
SPACE  ADMINISTRATION 
GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT  MD  20771 


ATTN  A  BRENNER  1 

ST  SYSTEMS  INC 
9701  L  PHILADELPHIA  COURT 
LANHAM  MD  20706 


ATTN  J  MAJOR  1 

ST  SYSTEMS  INC 
4400  FORBES  BLVD 
LANHAM  MD  20706-9008 


(1) 


NAVSWC  TR  90-409 


DISTRIBUTION  (CONTINUED) 


Copies 

ATTN  C  HANSON  1 

WORLD  DATA  CENTER  A 
FOR  GLACIOLOGY 

NATIONAL  SNOW  AND  ICE  DATA  CENTER 
CIRES  CAMPUS  BOX  449 
UNIVERSITY  OF  COLORADO 
BOULDER  CO  80309 


ATTN  A  OH  MURA  I 

GEOGRAPHISCHES INSTITUT  ETH 
WINTERTHURERSTRASSE  190 
CH  8057 

ZURICH  SWITZERLAND 


DEFENSE  TECHNICAL 
INFORMATION  CENTER 
CAMERON  STATION 

ALEXANDRIA  VA  22304-6145  12 


ATT X  GIFT  AND  EXCHANGE  DIV  4 
LIBRARY  OF  CONGRESS 
WASHINGTON  DC  20540 


INTERNAL  DISTRIBUTION. 


E211  iGREEN)  1 

E31  (GIDEPi  1 

E231  3 

E232  2 

K10  1 

K104  I DIDONATO)  I 

K104  (FELL)  15 

K105  2 

K106  1 

K 1 07  4 

K12  2 

K12  (CUNNINGHAM)  15 

K12  (GROEGER)  1 

K12  (SEAY)  1 

K13  2 

K14  2 

K40  1 

K41  (RUFTY)  1 

K43  2 

K43  (BROWN)  1 

K43  (SITZMAN)  5 


(2) 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMB  No  0704-0188 


Pubiu  reporting  burden  to*  this  collection  of  mtoi  motion  is  estimated  lo  average  1  hour  pei  response  indudinq  the  time  tor  reviewing  instructions.  seau hmg  existing  data  souru's 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  ol  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  ot  this 
collection  ot  information,  including  suggestions  tor  reducing  this  burden,  to  Washington  Headquarters  Services  Directorate  tor  information  Operations  and  Reports  12  IS  Jefferson 

Davis  Highway  Suite  1204  Arlington  V  A  22202  4  J02.  and  to  the  Office  ot  Management  and  Budget  Paperwork  Reduction  Project  10704  0188).  Washington.  DC  2050  f 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

July  1991 

3.  REPORT  TYPE  AND  DATES  COVERED 

Final 

4  TITLE  AND  SUBTITLE 

Forecasting  Geopotential  Trends  in  Southern  Greenland 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

James  P.  Cunningham  Patrick  J.  Fell 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Naval  Surface  Warfare  Center  (K12) 

Dahlgren,  VA  22448-5000 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NAVSWC  TR  91-409 

9  SPONSORING /MONITORING  AGENCY  NAME(S)  AND 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11  SUPPLEMENTARY  NOTES 

12a  DISTRIBUTION/ AVAILABILITY 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

1 3.  ABSTRACT  (Maximum  200  words) 

Various  data  types  including  SEASAT  and  GEOSAT  satellite  radar  altimetry  and  aircraft-collected  radio 
echo  soundings  were  used  to  develop  a  topographic  data  base  for  Greenland  south  of  72°  N.  The  data  base 
consists  of  the  ice  surface  and  underlying  bedrock  surface  elevations  in  regular  grids.  These  grids  and  given 
density  structures  provide  models  of  the  localized  mass  distribution  above  mean  sea  level  and  intermediate 
wavelengths  of  the  geopotential.  Contours  of  gravity  disturbance  and  deflections  of  t  he  vertical  at  altitude 
predicted  using  the  model  demonstrate  provinces  distinguished  by  their  gravimetric  gradients.  Comparisons 
between  modeled  gravity  disturbance  and  observed  free-air  anomalies  show  that  intermediate  wavelength 
trends  in  the  anomalies  are  often  predicted  by  the  model.  It  is  proposed  that  the  geopotential  model  could  aid 
in  the  design  of  an  airborne  geophysical  survey  of  southern  Greenland  and  in  the  interpretation  of  the  survey 
data. 

14  SUBJECT  TERMS 

SEASAT  GEOSAT  satellite  geophysical  data 

Greenland  gravity  topography 

IS.  NUMBER  OF  PAGES 

45 

16.  PRICE  CODE 

17  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

Standard  Form  298  (Rev  2-89) 


Orf,i,h.-.lh  iNSi  Strt  219  IS 


NSN  7540-01-280-5500 


GENERAL  INSTRUCTIONS  FOR  COMPLETING  SF  298 


The  Report  Documentation  Page  (RDP)  is  used  in  announcing  and  cataloging  reports.  It  is  important  that 
this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  its  title  page. 
Instructions  for  filling  in  each  block  of  the  form  follow.  It  is  important  to  stay  within  the  lines  to  meet 
optical  scanning  requirements 


Block  1.  Agency  Use  Only  (Leave  blank). 

Block  2.  Report  Date.  Full  publication  date  including 
day,  month,  and  year,  if  available  (e  g.  1  Jan  88).  Must 
cite  at  least  the  year. 

Block  3.  Type  of  Report  and  Dates  Covered.  State 
whether  report  is  interim,  final,  etc.  If  applicable,  enter 
inclusive  report  dates  (e  g  10  Jun  87  -  30  Jun  88). 

Block  4.  Title  and  Subtitle.  A  title  is  taken  from  the 
part  of  the  report  that  provides  the  most  meaningful 
and  complete  information.  When  a  report  is  prepared 
in  more  than  one  volume,  repeat  the  primary  title,  add 
volume  number,  and  include  subtitle  for  the  specific 
volume.  On  classified  documents  enter  the  title 
classification  in  parentheses. 

Block  5.  Funding  Numbers.  To  include  contract  and 
grant  numbers;  may  include  program  element 
number(s),  project  number(s),  task  number(s),  and 
work  unit  number(s).  Use  the  following  labels: 


c 

-  Contract 

PR  - 

Project 

G 

-  Grant 

TA  - 

Task 

PE 

-  Program 

WU  - 

Work  Unit 

Element 

Accession  No 

BLOCK  6.  Author(s)  Name(s)  of  person(s)  responsible 
for  writing  the  report,  performing  the  research,  or 
credited  with  the  content  of  the  report.  If  editor  or 
compiler,  this  should  followthe  name(s). 

Block  7.  Performing  Organization  Name(s)  and 
address(es).  Self-explanatory. 

Block  8.  Performing  Organization  Report  Number 
Enter  the  unique  alphanumeric  report  number(s) 
assigned  by  the  organization  performing  the  report. 

Block  9.  Sponsorinq/Monitorinq  Agency  Name(s)  and 
Address(es)  Self-explanatory 

Block  10.  Sponsorinq/Monitorinq  Agency  Report 
Number.  (If  Known) 

Block  11.  Supplementary  Notes.  Enter  information  not 
included  elsewhere  such  as:  Prepared  in  cooperation 
with...;  Trans,  of...;  To  be  published  in...  When  a 
report  is  revised,  include  a  statement  whether  the  new 
report  supersedes  or  supplements  the  older  report. 


Block  12a.  Distribution 'Availability  Statement . 
Denotes  public  availability  or  limitations.  Cite  any 
availability  to  the  public.  Enter  additional  limitations 
or  special  markings  in  all  capitals  (e  g.  NOFORN,  REL, 
ITAR). 

DOD  -  See  DoDD  5230  24,  "Distribution 

Statements  on  Technical  Documents." 
DOE  -  See  authorities 
NASA  -  See  Handbook  NHB  2200.2 
NTIS  -  Leave  blank 


Block  12b.  Distribution  Code 


DOD  -  Leave  blank. 

DOE  -  Enter  DOE  distribution  categories  from 

the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports 

NASA  -  Leave  blank. 

NTIS  -  Leave  blank. 

Block  13.  Abstract.  Include  a  brief  (Maximum  200 
words)  factual  summary  of  the  most  significant 
information  contained  in  the  report. 

Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 

Block  15.  Number  of  Pages.  Enter  the  total  number 
of  pages. 

Block  16.  Price  Code.  Enter  appropriate  price  code 
(NTIS  only) 

Block  17.-19.  Security  Classifications.  Self- 
explanatory.  Enter  U.S.  Security  Classification  in 
accordance  with  U.S.  Security  Regulations  (i.e., 
UNCLASSIFIED).  If  form  contains  classified 
information,  stamp  classification  on  the  top  and 
bottom  of  this  page. 

Block  20.  Limitation  of  Abstract.  This  block  must  be 
completed  to  assign  a  limitation  to  the  abstract. 
Enter  either  UL  (unlimited  or  SAR  (same  as  report). 
An  entry  in  this  block  is  necessary  if  the  abstract  is  to 
be  limited.  If  blank,  the  abstract  is  assumed  to  be 
unlimited. 


Standard  Form  298  Back  (Rev  2-89) 


